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Abstract 

We study correlators of 7£-symmetry currents in the Coulomb branch of Af = 4 super- 
symmetric gauge theory in the large- N limit, using the AdS/CFT correspondence. In 
particular, we consider gauge fields in the presence of gravity and scalar fields parameter- 
izing the coset SL(6, TR)/SO(6) in the context of five-dimensional gauged supergravity. 
From a ten-dimensional point of view these backgrounds correspond to continuous D3- 
brane distributions. We find the surprising result that all 2-point functions of gauge 
currents fall into the same universality class, irrespectively of whether they correspond 
to broken or unbroken symmetries. We show that the problem of finding the spectrum 
can be mapped into an equivalent Schrodinger problem for supersymmetric quantum 
mechanics. The corresponding potential is the supersymmetric partner of the potential 
arising in studies of the spectrum for massless scalars and transverse graviton fluctuations 
in these backgrounds and the associated spectra are also identical. We discuss in detail 
two examples where these computations can be done explicitly as in the conformal case. 
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1 Introduction 



For several years the dynamics of branes in string theory have been a fruitful playground 
to test strong coupling physics of gauge theories. For instance, the AdS/CFT correspon- 
dence |2], [| provides us with precise prescriptions to calculate correlation functions, 
spectra of gauge invariant operators, Wilson loops and c-functions in M = 4 super- 
symmetric Yang-Mills (SYM) theory in four dimensions at large N and large 't Hooft 
coupling. The data obtained this way from supergravity can sometimes be compared 
with field theory or provide non-trivial predictions for strongly coupled field theories. 
This correspondence can be extended also to theories with spontaneously or manifestly 
broken superconformal symmetry. Such theories arise either by giving vacuum expec- 
tation values to fields 0, p|-|[L3|| or by deforming the conformal theory with relevant 
operators [14|-|29|. Many of these deformations can be treated efficiently in the context 
of five-dimensional gauged supergravity ^Oj, |3l| and the resulting backgrounds have four- 
dimensional Poincare invariance and approach AdS§ in the ultraviolet (in a field theory 
terminology). Typically, towards the infrared, singularities appear which are not fully 
understood and seem to require a proper inclusion of the string theory dynamics or the 
use of other methods developed in gravity. 

In this letter we study correlation functions of 7?.-symmetry currents using the holo- 
graphic description of large- N gauge theories. For the conformal case correlation func- 
tions for operators in various representations of the ^-symmetry group S77(4) ~ £0(6) 
have been worked out in great detail (see, for instance, [^, |33|). Less is known about 
correlators in deformed gauge theories which are described by more general domain wall 
solutions of gauged supergravity. So far mainly scalars have been studied, namely the 
minimally coupled scalar ||, [K], [34[] (which has the same equation as the transverse trace- 



less graviton modes ||35||), active and inert scalars which parameterize deformations of the 
S 5 [^, |37], |38| , but also fermionic and abelian vector field fluctuations for the M = 1 
flow of |20| and the M = 4 Coulomb branch background of |J [lOj have been considered 



recently in |38 



We will show that for a specific class of examples this analysis can be extended to 
include fluctuations of non-abelian gauge fields which are dual to 7?.-symmetry currents 
of the gauge theory. We make a general connection between the fluctuation equation and 
supersymmetric quantum mechanics and find that, the relevant Schrodinger potential, 
associated with the spectrum, is just the supersymmetric partner of the potential arising 
from the corresponding massless scalar and transverse graviton-fluctuations equations. 
We show also that the corresponding spectra are identical. It seems plausible to us that 
this can be extended to the full set of fields in the supergravity multiplet. Using the 
AdS / CFT correspondence we calculate two-point functions of the symmetry currents in 
M = 4 SYM on the Coulomb branch in two particular cases.Q As expected, we find 
deviations from the conformal 1/r 6 fall-off for large separations r. From the non-analytic 

1 Other studies of the Coulomb branch of the AT — 4 SYM theory using the AdS/CFT correspondence 
can be found in 139(1 . 
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part of the correlator in momentum space we get contributions that are suppressed ex- 
ponentially for large separation. 

The choice of a particular state on the Coulomb branch breaks the 7?.-symmetry to 
a subgroup and therefore one might expect that broken and unbroken currents behave 
differently and in particular one would expect Goldstone bosons corresponding to the 
broken symmetry. From the dual supergravity point of view this symmetry is a local 
gauge symmetry and the massless bosons simply get eaten by the gauge fields and make 
them massive via the Higgs mechanism. Although the equations for broken and unbroken 
currents look quite different — they correspond, respectively, to massless gauge fields in 
a curved background and massive gauge fields — the associated spectra are identical. 
This result is not too surprising since on the Coulomb branch only conformal symmetry 
is broken but the currents still reside in the same supersymmetry multiplet. However, 
a small puzzle remains since the correlator has also an analytic piece that depends on 
which of the broken or unbroken currents are considered. For the two-point function of 
scalars such analytic terms give rise to contact terms and are usually dropped, but in the 
case of gauge field correlators they give rise to terms of the form x^Xv/r 6 , which might 
be interpreted in field theory as arising from Goldstone bosons. However, we do not find 
a one to one relation between broken currents and the presence of these terms in the 
correlators. We believe that these analytic terms are unphysical, since the corresponding 
mode is non-normalizable, and should be dropped. 

The organization of this paper is as follows: In section 2 we present some background 
material on gauged supergravity and calculation of correlators in AdS/CFT. We also 
make a general connection between the fluctuation equation and supersymmetric quan- 
tum mechanics. In section 3 we focus on our two main examples where calculations 
can be performed explicitly. We obtain the exact fluctuation spectrum of gauge fields, 
and the two-point functions in momentum and position space. In section 4 we give a 
summary of our results and give some final remarks. 



2 Generalities 

Our starting point is a specific truncation of the Af = 8 gauged supergravity action 



3lj including 5*0(6) gauge fields antisymmetric in i,j, with field strength F~~, 
where ju, v — 1,2,3,4,2; unhatted indices — 1,2,3,4 will be used later to denote 
Euclidean directions along the boundary at z — 0. For notational convenience we will 
occasionally use the collective index a = 1, 2, . . . , 15 to denote the adjoint representation 
of SO (6), instead of % and j or we will omit such an index all together. Furthermore, 
scalars in the 20' are represented by a symmetric traceless matrix M y . The action of the 
supergravity truncated to these fields has been constructed in j|(J and we follow closely 
their conventions. 

The Lagrangian density for the relevant fields of five- dimensional gauged supergravity 
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IS 

-^scalar ^ gauge > (1) 

where £ sca iar refers to the pure gravity-scalar sector and £ gauge contains the gauge fields 
and their interaction with the scalars and gravity. We first recall some results for the 
pure gravity-scalar sector since we are interested to study fluctuations of the gauge fields 
in the background of specific solutions of the gravity-scalar sector. The explicit form of 
the Lagrangian is 

^scalar = \^~ ^ {d^M d» M M' 1 ) - P , (2) 

where the potential is 

P = - 9 - [(TrM) 2 - 2Tr(M 2 )] , (3) 

with g being a mass scale. Alternatively we may use the length scale R via the relation 
9 = 2/R. 

Super symmetric solutions of (||) preserving 16 supercharges and Poincare symmetry 
in four-dimensions have been studied extensively and they correspond to states on the 
Coulomb branch of M = 4 SYM theory. Their interpretation in ten dimensions is simply 
in terms of a continuous distribution of D3-branes. For these backgrounds the matrix of 
scalar fields can be brought to a diagonal form using a gauge transformation. Thus we 
are left with six scalar fields that parameterize 

M = diag(e 2/ \...,e 2/36 ) , (4) 

obeying the constraint Y2^=i A = 0- There are five independent scalar fields, denoted by 
ai, I = 1, 2, ... 5, and the relation to the fa's is given by fa = Yli=i ^u a ii where Xu is 
a 6 x 5 matrix, with rows corresponding to the fundamental representation of SL(6, H); 
the normalization conventions can be found in eq. (2.4) of [03]. The metric ansatz reads 



ds 2 = e 2A{z) {dz 2 + rj^dx^dar) = dr 2 + e 2A{r) r]^dx^dx u , (5) 

where the relation between the coordinates z and r is such that dr = —e A dz. In addition, 
all scalar fields depend on the variable r or equivalently z. The most general solution 
preserving 16 supercharges has been found in and is conveniently presented in terms 



of an auxiliary function F(g 2 z). Specifically, the conformal factor is given by 



e 



2A 



9 2 {-F'f" , (6) 



where the prime denotes the derivative with respect to the argument of F(g 2 z). In 
addition, the profiles of the scalar fields are 



em = J~br / = rK F - 6 *)> i=l,2,...,6. (7) 

i=i 
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The constants of integration are ordered as b\ > 62 > • • • > b e and the function F is 
constrained to obey the differential equation 

{F'Y = f. (8) 

Equating n of the integration constants hi (or equivalently the associated scalar fields /%) 
corresponds to preserving an SO (n) subgroup of the original 50(6) 7?.-symmetry group. 
We note in passing, that there is a deep connection between solutions of the gravity-scalar 
sector of the five-dimensional gauged supergravity that we just reviewed, and the theory 
of algebraic curves and associated Riemann surfaces to which the differential equation 
(D is related |1J, |TJ. 

Let us now turn to the part of the Lagrangian containing the gauge fields. First, we 
have to replace the partial derivatives in @ by gauge- covariant ones d^M 1 ^ — > ^M* 3 + 
g(A~M kj + A~ t M lk ), and, second, we add the gauge kinetic term 
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£ gauge = --{M-'fiM-^F&F^ . (9) 



Since we are interested in two-point functions we only need to keep terms in ([!]) and 
(P) which are quadratic in the gauge fields and the scalar fluctuations in the symmetric 
unimodular matrix M. Note that although for our solution the matrix M is diagonal as 
in (|j), we have to consider fluctuations along the diagonal as well as off-diagonal ones. 
Using the fact that M is diagonal (§) for our backgrounds, we collect all terms that can 
give quadratic terms in the fluctuations of the scalars and the gauge fields 

^Auad. = -^-^^^^-^sinh 2 ^-^)^^^ 



--Tr ((d^MM- 1 - M-%M)A^) (10) 

8 V J quad. 



-^Tr^MM-^MM- 1 ) - P| quad . . 

The first line above is already quadratic in the gauge field fluctuations. We emphasize 
that F~L = dfiA- —d„A~ is, for our purposes, the relevant part of the gauge field strength. 
The second line in the above expression is already linear in the gauge field fluctuaction. 
Hence, we are supposed to expand it to linear order in the scalar fluctuations. Finally, 
the third line has to be expanded to quadratic order in the scalar field fluctuations. In 
this paper we are only interested in the gauge field fluctuations which, however, couple 
to fluctuations of the scalars. Therefore, it is not a priori correct to simply keep the 



terms in the first line in ( jTOj ) and drop the rest. Nevertheless, we will now explain that 
this procedure gives the correct result since there is a field redefinition that effectively 
decouples the gauge field fluctuations from those of the scalars.0 To see that let us expand 

2 We thank M. Bianchi for prompting us to explain in detail how the decoupling between scalar and 
gauge field fluctuations actually works as well as for other related comments. 
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the second line in ( fL0|) and keep the linear term in the scalar field fluctuations. We find 
that 

-§Tr Ud n MM' 1 - M' l d n M)A^\ 

o V / quad. 

= |((e- 2ft - e- 2 ^)d^M l3 + 2(e- 2ft «9^ - e~ 2 ^ d^ t )bM^ A% . (11) 

From this we immediately deduce that the diagonal fluctuations SM^ do not couple to 
the gauge fields. A less trivial fact is that the scalar fluctuations in 5Mij that belong to 
any unbroken subgroup of 5*0(6) do not couple to the gauge fields as well. The reason 
is that in this case $ = @j, since then the corresponding integration constants in (|7|) are 
equal, i.e. 6j = bj. Hence, let us consider the remaining cases with /3, ^ (3j which arise 
when the indices i,j belong to the coset. If we make the field redefinition 

4'-4+^(^). a* ft. (12) 

the mixed terms between scalar and gauge field fluctuations in (p~Q|) (with the substitution 
(|TT|) understood) disappear and the fluctuations decouple. Note that the field redefinition 
(p"2|) acts as an abelian gauge transformation and as such it leaves the gauge field strength 
invariant (to the quadratic order we are working). We emphasize that the field 
redefinition (|T^) does not guarantee that there will be no mixing between scalar and 
gauge field fluctuations at the cubic or at some higher order in the fluctuating fields, but 
only that the quadratic fluctuations decouple. We also note that a similar decoupling 
mechanism for vector and scalar fluctuations was found to be at work for the flow of |^D[ 
in pg| . There, it was observed that decoupling was achieved since the gauge field and a 
(charged) scalar appeared in a gauge invariant combination. 

The field redefinition (fl2|) removes the scalar fluctuations of 5Mij since it removes 
terms quadratic in first derivatives of SMij from the Lagrangian. The remaining terms 
are at most linear in first derivatives and of the form B^bM^bMij + B^-bM^d-^bM^ for 
some space-depended B^ and B^ which are symmetric in Clearly the derivative- 
term can be removed by adding an appropriate total derivative so that we are left with 
a non-dynamical field SM^ corresponding to no physical degrees of freedom. What we 
have is nothing but a manifestation of the Higgs effect in a curved background. As in 
flat space-time, the Goldstone bosons corresponding to the broken gauge symmetries are 
eaten by the gauge bosons which then become massive. 

Since we are only interested in the gauge field fluctuations we ignore the scalar fluc- 
tuations for the rest of the paper and concentrate on those for the gauge fields which, 
after the redefinition (|12l), are described by the first line in (|I0|) 

^C(A) quad . = -\^ 2{ ^F%Fff - £ sinh 2 (A - &)4'4 • (13) 

The second term corresponds to mass terms for the gauge fields, if the scalar fields f3i are 
not equal. This implies that for general states on the Coulomb branch the bulk gauge 
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symmetry S0(6) is spontaneously broken and, hence, that the 7?.-symmetry group of the 
field theory on the boundary is reduced accordingly. Notice also that the kinetic term 
for the gauge fields is not canonically normalized as it gets "dressed" by the scalar fields. 
This will have important consequences, as we will see. 

The equation of motion following from this quadratic action ( |T5| ) is: 

bA% : D n (e- 2 ^+^F>? ) - g 2 sinh 2 (A - ^)A% = . (14) 

In solving these equations we have to distinguish two cases: First, for the unbroken 
symmetry (currents), for which $ = /3j, we can use the gauge symmetry to choose the 
gauge A % i = 0. This still allows for restricted gauge transformations with parameters 
that depend only on the x M 's, but not on z. Then, the p, = z component of the eqs. 



(14) yields the constraint d z d^A^ = which allows to eliminate unphysical longitudinal 
modes via a restricted gauge transformation. The equation of motion for the remaining 
physical (transverse) modes Aj^ which obey d^A^ = is the same for all components 
and can be written as an equation for a scalar field, which we denote by $: 

d z (e B d z $) +m 2 e B $ = , (15) 

with the definition 

B = A - 2(J3i + fa) . (16) 

To arrive at this equation we have performed a Fourier transform in the a^-directions 
with = — m 2 .[] 

For the broken symmetry currents for which $ ^ j3j we cannot use a gauge symmetry 
to eliminate degrees of freedom. In order to calculate the two-point functions we couple 
the gauge field to an external source by adding —\A l yJ^ to the gauge field action flT3|). 
The source is required to be covariantly conserved, i.e., D^J~ = 0. We choose to decom- 
pose the gauge field into transverse modes Ah, longitudinal modes <9 M £ = A^ — Ah, and 
the component A z . The equations of motion (|I4"D give 

d z (e B (d z Ai-d ll A z + d z d,0) + e B nA^ - e°{UA^ + d^) = e 3A J z (17) 

and 



e B (DA z -d z nZ)-e c A z = e 3A J z , 



where □ = r]^d^d u . The above coupled system of equations can be further simplified. 
By taking the derivatives and d z in (|T7|) and ( fHf) respectively, adding up the resulting 
expressions and then using the condition D^J^ = 0, we obtain a relation that determines 
£ in terms of the component A z , namely 

e c n£ + d z (e c A z ) = , (19) 



3 For notational simplicity we did not include indices i,j in defining B in (|l6|). Nevertheless it should 
be kept in mind that different choices for the scalar fields (3i and j3j lead to different values for B. 
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where 

e c = g 2 e 3A sinh 2 (& - fi s ) = ~g 2 (k - ft^V* . (20) 

The first equality defines C, whereas the second one follows with the help of (|7|) and relates 
C to B which was defined in QTEp. Using ( [TUj ) to solve for □£ and then substituting back 
the result into fll8f ) we find the equation for the mode A z , which decouples from the 
transverse modes: 

e B DA z + e B d z (e~ c d z {e c A z )) - e c A z = e 3A J z . (21) 



With further manipulations using (21), we may cast (|17]) into an equation for the trans- 
verse modes 

e B DAi + d z {e B d z Ai) - e c 'A$ = e* A J± , (22) 

where we have defined the transverse current-source as = (5^ — <9 M <9„/D) J^. In 
order to compute the two point functions in momentum space we need solutions of the 



homogeneous equations ([H]) and (|22f) . Actually, after a Fourier transform in the x M 



brane-directions, we can write both equations as an equation for a scalar field 

d z (e B d z $) + [m 2 e B - \g\b t - b^A $ = , (23) 

where we have dropped the source term. Its effect will be implemented by imposing 
appropriate boundary conditions to the solutions. For the case of Q22D the scalar $ 
denotes any component of Ah. In order to cast ([2~Tf) into the form (p3|) , we have used 
(p0|) and defined $ = e c A z . For /3j = f3j we recover from (^3|) eq. (|15|) that describes the 



cases with unbroken symmetry. Hence, for full generality, we may use (|23|) in order to 
calculate current- current correlators. We will follow the standard procedure of P, 0] and 
we will work in Euclidean signature unless stated otherwise. 



In order to proceed we need a complete set of eigenfunctions of (23), which for the 
examples we will discuss in the next section can be found explicitly and is given in terms 
of hypergeometric functions. Furthermore, we keep the solutions that blow up at the 
AdS boundary since they correspond to current operator insertions || ||. Finally, we 
have to evaluate the on shell-value of the action -\ f d 5 x£ with -K = for solutions 
$ of (|23|).^ We find the boundary term 



N 2 

- lim -e $<9 2 $ 

e-0 327T 2 



N 



2 



167T : 



k 2 H(k) . (24) 



In order to keep formulas short in later sections we have written out the overall factor 
1/k 2 in the definition of H(k). In order to obtain the correct result we have to normalize 
*&\z=e = 1 and take the limit in fl24]). Re-introducing Lorentz and group theory indices 



4 The overall normalization is found by carefully keeping track of all the prefactors in the dimensional 
reduction in the 5' 5 -directions of the ten-dimensional type-IIB action to five dimensions. In particular, 
73 = j^R S - Tnen using 2n 2 w = (2n) 7 a^gl, R A = Ang s a' 2 N and ^55 = tt 3 wc find the result mentioned 
above. 
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properly, we can present the current-current correlators in momentum space schematically 

as 

N 2 f k k \ 

(J a »(k)4(-k)) = ^- 2 S ab (V - -j^J k'G{k) , (25) 

where a group theory factor and the momentum space version of the projector, which 
guarantees that the amplitude is transverse, have been included. The factor H(k) = 
k 2 G{k) depends also on the adjoint indices a, b, but for reasons similar to those explained 
in footnote 2 we have not explicitly displayed them. 

In the explicit calculations performed later in section 3 we will not use H(k) directly, as 
defined in (0), because the correlator in x-space is too singular to be Fourier transformed 
to momentum space. However, by using differential regular izat ion one can make sense of 
such expressions by writing singular functions as derivatives of less singular ones and then 



defining the Fourier transform by formal partial integrations [fjlJl . In our case we have 
to take the correlator to be of the form ~ □□G(x) which is just k 4 G(k) in momentum 
space. Hence, the correlator in x-space becomes 

N 2 

{J a ^)J») = ^- 4 5 a \nS, u - d,d u )UG{x) , (26) 



where 



i\i,Ak-x H \k) _ 1 



oc 



G(x) = —r / dV- x — = - / dkH(k)J x (kr) , (27) 
vk 1 I k 







with Ji(kr) being a Bessel function. 



2.1 Supersymmetric quantum mechanics 



In this subsection we want to study general aspects of the fluctuation equation (|23|), 
before we proceed in section 3 to describe two special cases where calculations can be 
performed exactly. Writing $ = e~ B//2 \l/ the field equation (p3|) turns into the one- 
dimensional Schrodinger equation 



m 2 * 



(2f 



with potential 



V 



im'f + -B" 
4 V 1 2 



(29) 



This potential, though not at all obvious, can be cast into a form that appears in super- 
symmetric quantum mechanics. First, we rewrite it differently using the properties of 
our solution (H) and (^) and in particular (H) which proves useful in turning derivatives 
with respect to the variable z into functions of the auxiliary function F only: 



V 



4 fl/2 



g'f 



64 



•E 



1 



E 



i=i 



(30) 
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Comparing with eq. (4.16) of JL3| (after setting in there the parameter A = 4) we find 
that this can be written solely in terms of the conformal factor in the metric ansatz (^]) 

V = ? a' 2 - - A" . (31) 
4 2 y ' 

This potential^ has the same form as the potential appearing in supersymmetric quantum 
mechanics |42|, f|3[ with superpotential W = —3/2 A'. In fact, it is the supersymmetric 
partner of the potential 

V. = \a' 2 + \a» , (32) 

that appeared in studies of 2-point functions for scalar fields or transverse graviton fluc- 
tuations || 35, 11, 12, 44, the relation of ( |32l) to supersymmetric quantum mechanics 
in the context of gauged supergravity was first hinted in [11 and explicitly noted in [44 



Note that, the Schrodinger problem is universal and does not depend on the indices 
i,j of the gauge currents. Consequently, the mass spectrum is the same irrespectively of 
whether it is associated to currents corresponding to broken or unbroken symmetries. In- 
stead, the wavefunction $ does depend on the indices i,j through the explicit dependence 
on them of the conformal factor B defined in ([TB]) (cf. footnote 3). 

It is well known from the general theory of supersymmetric quantum mechanics that 
the spectra of superpartner potentials, such as fl3~T|) and (|32|), are identical except for a 
zero mode. However, in our case such a mode is not normalizable due to the asymptotic 
behavior of the function A(z) as z — > and, therefore, is not included in the spectrum. 
Hence, the spectra of current fluctuations, corresponding to ([H]) and those for dilaton and 
transverse graviton fluctuations, corresponding to (|32|), exactly coincide, as advertised in 
the introduction. We note, that related observations concerning a 50(3) invariant sector 



of 5d gauged supergravity and a particular Coulomb branch flow have been made in [[38 

B 

The analysis of the qualitative features of the spectrum can be done in a similar 
fashion as in the case of the superpartner potential arising in the case of scalar correlators 



fT|, [12| . At the boundary z = the potential goes to +oo as V — A. The behavior 



in the interior depends on the number n of constants of integration bi that equal the 
maximum constant among them, b\. We follow closely the discussion of [|ll|, |13| to which 
we refer for further details. For n = 4, 5 the range of z necessarily extends to +oo, i.e. 
< z < oo, corresponding to F = b\. We find that, for n = 5, the potential goes to zero 
as z — > oo and the spectrum is continuous. For n = 4 the potential approaches a constant 
value, as z — > oo, which is given by V m i n = ^r/ Q • Therefore, although the spectrum 
is continuous, there is a mass gap whose squared value is given by the minimum of the 
potential. For n = 5 the potential behaves as 

K T/ 15 / 4 

n = 5 : V 5 ~ — — , as z —>■ oo . (33) 
z l 



5 An alternative way to prove the equivalence of the potentials ( |29| ) and (|3j) is to use the differential 
equation obeyed by the /Vs, namely (i[ = A' + |e A+2ft 

6 The authors of p8| informed us that their arguments concerning graviphotons are actually broader 
and include all massive cases where U(1)r is broken. 
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For n = 1,2,3 the potential goes to +00 as F — > b\ and therefore the spectrum must 
be discrete. Therefore there should be a maximum value for z, denoted by z max , that is 
determined by solving the algebraic equation F(z max g 2 ) = b\. We find the behaviour 



n 



1,2,3 



Cn 



(z - Z ) 2 



(4 - nf 



1 

4 



as z 



(34) 



For more details on the full structure of the potentials ( pi]) and (|32|), which generically 
can be written using elliptic functions, the reader is referred to the original literature 
In the two special cases, to which we turn now in section 3, all computations 



11, is 



and results can be written in terms of elementary functions. 



3 The 2-point function 

In the previous section we introduced all necessary ingredients for the calculation of 
correlators of symmetry currents and pointed out the relation between supersymmetric 
quantum mechanics and the fluctuation equations. In this section we want to use these 
results and apply them to two specific backgrounds worked out in H, |T(], [TTJ. These 
backgrounds correspond to distributions of D3-branes on a disc or a three-sphere 0, Q 
and they both break the bulk gauge symmetry down to 5*0(2) x 5*0(4). The broken 
symmetries form the coset S o{i)Tso(4) • ^ n ^ e dual ne ^ theory side these backgrounds 
correspond to states on the Coulomb branch of M = 4 SYM theory with reduced 1Z- 
symmetry. In the following we will calculate the correlators in momentum and position 
spaces. 



3.1 Distribution of D3-branes on a three-sphere 

We begin our exactly solvable examples with the case of a model representing D3-branes 
uniformly distributed on a three-sphere. The expressions for the metric and the scalar 



fields have been given in || |TT|] . The five-dimensional metric (|5|) has the conformal factor 

r„ COS 2 / 3 W n 7T 

R 2 sm 2 u 2 



e 2A = -^— 2 — , < u < - , (35) 



where we have defined for notational purposes the dimensionless variable u = r^z/R 2 . 
The parameter r actually plays the role of the radius of the three-sphere. The AdS 5 
boundary corresponds to u — 0, whereas at u = it/ 2 there is a naked curvature singularity. 
This is however naturally interpreted, from a string theoretical point of view, as the 
location of the distribution of the D3-branes on the three-sphere. 
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The profiles of the scalar fields are 

e 2Pi = e 2/3 2 = cog -4/3 u ^ e 2/3 3 = = £ 20 e = ^2/3 u ^ 

From a ten-dimensional view point, these scalars deform the five-sphere line element that 
appears in the D3-brane solution in such a way that the subgroup 50 (2) x 5*0(4) of the 
isometry group 50(6) is preserved. The Schrodinger potential (|TL|) is found to be 

R \ sin 2u J 

It is not difficult to show that a complete orthonormal set of solutions to the corresponding 
Schrodinger equation is given by 



%i= l2Mn^^cj^u , „=i, 2 ,..., (38) 

V n + 1 R 2 sm L/2 u 2 

where the Pn 's are Jacobi polynomials, provided that the spectrum is given by 

4r 2 

m 2 n = ^n(n + l) , n = l,2,... , (39) 

Note that the case with n — 0, giving rise to a zero-mass eigenvalue, is not included in 
the spectrum since the corresponding Schrodinger norm diverges. The eigenvalues 



coincide with those found for dilaton fluctuations in || [HJ using the same background as 
here, in agreement with our general discussion in section 2. Also the n-dependent overall 
constant in (|38|) has been chosen such that the \l/ n 's are normalized to one. 

The conformal factor appearing in the equation of the fluctuations (|23|) is: 



cos° u 
sin u ' 



R 



sin u cos u 
cos u 
sin u ' 



*,J = 1,2 
= 3,4, 

1,2, j = 3,4,5,6. 



e B = — x I 1 i 7 = 3 4 5 6 (40) 

D I sin?/, cos ii. V / 



3.1.1 The 2-point functions 



Using (^), ( |36D and (|40"D we find the wave equation for the transverse modes of the gauge 
field in the unbroken 50(2) subgroup, the coset and the unbroken 50(4): 

(i-x)(x 2 &y- ^$ = o, 

(I - x)(x&)' - - = 0, x = cos 2 u G [0, 1] , (41) 

12 

x(l - x)$" - — $ = 0, 

where the prime denotes derivatives with respect to x and k 2 = R 4 /r^k^k^ 1 , i.e. is the 
length-square of the four- vector k^ rescaled for notational convenience with the indicated 
factor. 
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The wave-functions that blow up at the boundary at x — 1 and are regular at the 
singularity at x = are given in terms of a hypergeometric function as[] 

$ = T((3 + A)/2)r((3 - A)/2)x x F ^A, 2 , x] . (42) 



where A = \/l — A; 2 , and where we have introduced the parameter A = 0, | and 1 for the 
currents corresponding to the unbroken SO(2), the broken coset and the unbroken 50(4) 
symmetries, respectively. The proportionality constant in (f4~2"|) has been fixed such that 
$(1) = 1 and hence at the boundary the solution becomes proportional to a ^-function, 
i.e., fully localized operator insertion. It is interesting to note that the wavefunctions $ 
in all three cases differ only by different powers of x. This is related, as we have seen, 
to the fact that the mass spectra for broken and unbroken currents are identical. From 
(f42|) we extract 

H(k) = + l/4(*((l + A)/2) + *.((l-A)/2) + 2 7 ) 

= A If, 2n + P 

k 2 2 ^ n(4n(n + 1) + k 2 ) 

which has a discrete spectrum of poles at k 2 = —4n(n + 1), n — 1, 2, . . ., corresponding 
precisely to the mass eigenvalues fl3~9|) . However, if A ^ 0, there is an additional pole at 
k 2 = 0. We will comment on this in various places below. 

The three correlators differ only in the coefficient of the 1/k 2 term. In the case of 
scalar correlators this would just give a contact term and could be ignored, but in the case 
of the symmetry-current correlators this has important consequences as we will explain 
shortly. Using fl2T|) we obtain the following exact expression for the function G(x) in the 
correlator (p6|): 

where K± denotes the modified Bessel function and in writing the term containing lnr 
we discarded an infinite constant. We have also dropped a 1/r 2 term, which, since 
□1/r 2 ~ 5^'(r), contributes only contact terms to the correlator which we consistently 
ignore. Hence, we find 



2 o 3 00 

DG(x) = A^ + ^^(2n + l)^(^^i(2^Tl)^) • (45) 

n=l 



Let us perform the consistency check that for small r, or equivalently, in the limit 
r — > 0, we should recover the conformal result. The dominant contribution in this limit 



7 Throughout the paper we will make use of special functions and their properties following the 
conventions of lEEl. 



12 



comes from the infinite sum which can be approximated by an integral 

1 f 1/r dn 1 
G(x) = — — + . . . ~ -— lnr 2 , as r^O. (46) 
2r z Jx n 4r z 

This gives rise to 

DG{x) ~ , as r -> , (47) 

which in turn, gives a 1/r 6 fall off for the correlator ( P6|) at short distances. As expected, 
this coincides with the result in the conformal case (see, for instance, eq. (30) of f32|). 

The behavior of G(x) for large r is easily found from the asymptotic expansion of 
the modified Bessel function. For large r each separate term in the infinite sum behaves 
as e - m n r / / r 3 / 2 ) where m n are the mass eigenvalues in (|39|) and hence gives rise to an 
exponential fall off. Keeping the two most dominant contributions in the right hand side 
of (HD we obtain 

G(x) ~ A-^hir + !v| (&\ 3/2 e -2V2r 0r/R i as r ^ ^ . (48 ) 
2-R 4 8V2 V r o r / 

For the cases corresponding to the broken coset currents and the unbroken SO (4) currents 
we have A ^ and therefore the dominant contribution for large r comes from the first 
term in (lilf) . When substituted into the correlator in ( pop it produces a contact term, 
which we drop, and a term of the form 

N 2 r 2 1 

< w jj(o)> * X5ab ^it^ {rH ^ ~ Ax » Xv) ' as r ^ 00 • (49) 

This term decays only with the forth power of the distance and at first sight it might be 
tempting to interpret it as arising from the massless Goldstone boson associated with the 
broken symmetry^ From a physical point of view there are several problems with such 
an interpretation: First, this term does not appear on equal footing for all three types of 
currents although they reside in the same supersymmetry multiplet. Its existence might 
seemingly be acceptable or even desirable for the broken symmetry, but this term also 
appears for the unbroken SO (4)-symmetry currents. We also know from section 2 that 
the gauge fields dual to the broken currents become massive via the Higgs mechanism 
and, therefore, are not expected to produce any massless states. Second, the pole of 
the massless state corresponds to a non-normalizable mode and it is not expected to 
show up in the two-point function. The most plausible solution seems to be that these 
poles are actually unphysical and should be dropped from the correlators. Note that a 
similar problem was found in |36| for the two-point function of active scalars in the same 



backgrounds we are discussing here. The mysterious massless poles in that paper were 



later shown to be absent if a different prescription for the correlators is used f37fl . It 
seems likely, although we have not checked, that an improved prescription would resolve 
the puzzle in our case as well.[] 

8 Work on the AdS/CFT correspondence and the Goldstone bosons has been reported using a different 
model in ^jj . 

9 Actually, we were able to explain the presence of these massless poles we found in the supergravity 
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3.2 Distribution of D3-branes on a disc 



Our second exactly solvable model represents D3-branes uniformly distributed on a disc 
of radius r . The expressions for the metric and the scalar fields have been given in 
|| [11] . The five-dimensional metric (|5p has the conformal factor 



9 A ri cosh 2 ^ 3 u , , 

e 2A = -± . - 2 , 0< M <oo, (50) 
ti smh u 

where as before u = r z/R 2 . The scalar fields are given by 

e 2 ^ = . . . = e 2/34 = cosh 2 / 3 u , e 2/3s = e 2ft = cosh" 4 / 3 u . (51) 

As before, from a ten-dimensional type-IIB view point, these scalars deform the five- 
sphere line element that appears in the D3-brane solution in such a way that the subgroup 
SO (2) x 5*0(4) of the isometry group 5*0 (6) is preserved. The Schrodinger potential ( pT|) 
becomes 

V = ^ (l + ) . (52) 

R 4 V sinh 2 2u) V ' 

The energy spectrum for this potential is continuous and has a mass gap 

m 2 > -| ■ (53) 

As before the zero mode corresponds to a non-normalizable wavefunction. 
The conformal factor appearing in the equation of the fluctuations fl2"5| ) is: 

f , i,j = 1,2, 

sinh u ' J 

eB = i x \ sium. U = 3,4,5,6, (54) 
I TO. * = 1 ' 2 ' J = 3,4,5,6. 

3.2.1 The 2-point function 

The wave equation (p3|) for the gauge fields of the unbroken 5*0(2), the broken coset and 
the unbroken 50(4) symmetries, respectively, are: 

x 2 (l - xW' - — <£> = , 
4 



k 2 1 
— $ ^_ 

i I cosh^ u 

k 2 



x(l-x)(x&)' -— $- 7 (l-z)$ = 0, x = — — e [0, 1] , (55) 



1 -x)(x 2 <$>')' - '—<& = 0, (56) 



calculation by a field theory calculation in the free field approximation. These results are added as an 
addendum at the end of this paper, since they were found after publication of the original version of the 
paper. 
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where, as before, k 2 = k 2 R 4 /r 2 ,. The properly normalized solution that is also regular in 
the interior is 



= r((l + A)/2)r((3 + A)/2) (1+A)/ / A-l A + l \ 

r(i + A) " V 2 ' 2 J 



(57) 



where and similarly to before, the parameter A = 0, ~ and 1 for the currents 

corresponding to 50(2), to the coset and to SO (4), respectively. From this we obtain 

H{k) = ^ + i(V((l + A)/2)+ 7 ) 
A-l 1 r , e -« - 



fc 2 2 /„ 1 - e- 



dt — (5* 



and then 



! ^ If 30 , y e -vV^W 



= -(l-A)-^lnr + — / dy^- (59) 

1 ' 21 'R A 2r 2 J y sinh yy/yt + ry/R* 1 ^ 

Using this result it can be easily seen that the short distance behavior of the propagator 
is the same as in the conformal case and in particular ( fjBp is recovered. At large distances 
one finds that the two most dominant terms are 

G(x) i(l-A)^lnr + Z^J^ e -nW* 2 as r ^ ^ ? (60 ) 

where naturally the range of the Yukawa-term is set by the mass gap in (0). Hence, for 
the case where A ^ 1, corresponding to the cases of the broken coset and the unbroken 
5*0(2) symmetries, the first term dominates for large r giving a contribution to the 
correlator similar to (ff5D, but with A replaced by 1 — A. For similar reasons to those 
that we outlined for the case of the sphere-distribution of D3-branes after (|4"9"D, the 
interpretation of such a term as being related to the Goldstone bosons is problematic 
and we believe that they are unphysical. (However, see footnote 9.) 



4 Discussion 

In this letter we studied ^-symmetry current correlators in certain states on the Coulomb 
branch of Af = 4 SYM using the standard description of the AdS/CFT correspondence. 
The surprising result is that the spectra derived from the analytic structure of the correla- 
tors agree with spectra of other operators corresponding to dilaton and to the transverse 
graviton fluctuations. Furthermore, it turned out that the spectra are identical and do 
not depend on whether they are in the unbroken part of the left over global symmetry 
or reside in the coset, except for certain zero-mass poles which do depend on the sector. 
These poles give rise to a 1/r 4 fall off of the correlators at large distances, the behavior 
expected of massless scalars, but we did not find good physical reasons to identify them 
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with Goldstone bosons of the broken symmetry currents. We rather think that these poles 
are unphysical since they correspond to non-normalizable states and are inconsistent with 
the fact that the currents are all in the same supersymmetry multiplet.Q 

Rephrasing the fluctuation equations into a supersymmetric quantum mechanics prob- 
lem we found that they all fall into the same universality class and, furthermore, the 
Schrodinger potential are the supersymmetric partner potentials arising from the dilaton 
or from the transverse graviton fluctuations, which are identical. This indicates that all 
fluctuations in such backgrounds fall into the same class of supersymmetric quantum 
mechanics problems. 

To obtain a more complete picture including the Goldstone bosons one probably has 
to include additional modes that live on the D3-branes which create the singularity in the 
infrared. In our set up with a continuous distribution of branes this seems a formidable 
task, and as a starting point it seems more feasible to study simpler examples, e.g. two 
stacks of coinciding branes or a single test brane separated from a stack of branes, in 
which case one would readily know the additional modes and their respective couplings 
to the bulk fields. We leave these issues for future work. 

It will also be interesting to investigate current-correlators using solutions of D = 7 
and D = 4 gauged supergravity that are dual to the (2,0) theories in six dimensions and 
the three-dimensional theories with sixteen supercharges, respectively, on the Coulomb 
branch. For a class of such backgrounds corresponding to a scalar-gravity sector analo- 
gous the one used in the present paper the most general solution has been found and is 
very similar to that in |13| (see also |I2"|| ). The spectrum, of fluctuations corre- 



sponding to a massless scalar has been also exhaustively studied and in some cases the 



computations can be performed explicitly |T3L Similarly to the present paper, in these 



cases as well, it is quite plausible that the current-correlators and the associated spectra 
are related via supersymmetric quantum mechanics to those of the massless scalar. 
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Addendum 

The purpose of this addendum is to investigate the structure of the massless poles that 
appear in 2-point functions of broken symmetry currents in M = 4 SYM theory using 
purely field theoretical techniques and to compare the results with those obtained in 
section 3 using supergravity and the AdS/CFT correspondence. We had completed the 
essential part of this work around March of 2001. Parts of it are based on ideas developed 
around that time in collaboration with D. Freedman and K. Skenderis. 



General formulation 

We start with the case of unbroken 7£-symmetry where the vev's corresponding to the 
six scalars of the theory are turned off. The ^-symmetry currents are represented as 
bilinears in the scalar fields X % , i — 1, 2, . . . , 6 transforming in the adjoint of SU(N) 

J- = -^-T^{X%X j ) + fermions , (61) 
9ym 
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where T a are 6x6 matrices of SO (6). The scalars X\ being free fields, obey the following 
two-point function Q 

(X; q {x)Xi s {0)) = g 2 YM ^(5 gr 5 ps - ^ pq 5 rs )^ , p, q, r, s = 1, 2, . . . , N . (62) 

After performing the Wick contractions we compute the two-point function for the cur- 
rents 

(Jfa)J b M) ~ NH a \U5, p - d,d v )- , (63) 
where we have kept only the leading term in the 1/A^-expansion.Q This is indeed the 



correct result for the two point function which also agrees with the AdS/CFT result [[32 

In the case that the symmetry is broken by turning on non-zero scalar vev's, we 
replace X % by X l vev + 5X l , where the SX l have the same free field two-point function 
as in (0). Besides the bilinear term (|6lD the current contains now a term linear in 
fluctuating fields 

SJ l = -^T^riX^d^) , (64) 
9ym 

where we have introduced the vev's 

N 

x; ev = (X' 1 ) = dia g (xj, xl,...,x i N ), J2 x l = - ( 65 ) 

p=l 

At this point it is convenient to replace the adjoint 50(6) indices by a = [ij] and b = [kl]. 
Then, the matrix elements of the S0(6) generators become T^ n = 5 im 5j n — Sj m 5 in . The 
leading order correction to the conformal result (|63|) is 



(5J^(x)5J k J(0)) ~ J-fP^d^i , (66) 
9ym r 



where the group theoretical factor i?*,- jy takes the form 



N 



H^ kl = b lk A 3l - 5 jk A a - 5 u A jk + b 3l A lk , A l3 = ^ X^ . (67) 

P =i 



It is clear that, in the UV where the vev's can be neglected, the conformal result ( |63"D 
dominates, whereas in the IR the dominant term is fl6"6|). The symmetric tensor Aij is 
given in terms of the scalar vevs only and depends on their distribution. In the following 
we think of the vevs X % vev as defining points in R 6 . In most examples we use the 
fact that in the large A^ limit such a discrete distribution can be well approximated by a 

11 In our conventions the field theory action has an overall factor of \j g\ M . 

12 For finite N, the 1/iV-term in ( |6^ ) induces a shift which replaces the coefficient iV 2 by iV 2 — 1 
corresponding to the dimension of the SU (N) group. We also note that the contribution of the fermions 
only affects the result by an overall ./V-independent numerical constant which is not important for our 
purposes. 
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continuous one.0 Furthermore, we will consider situations where the distribution spans 
only a lower dimensional submanifold embedded in R 6 . The tensor H^' kl contains all 
the important information about the zero mass poles. It is antisymmetric in the indices 
ij and kl separately and symmetric under pairwise exchange. Note that only if both 
indices i,j are along the vev-distribution Aij is non-zero. That implies that if all indices 
correspond to directions which are perpendicular to the distribution, then H^ ,M = 0. 



Basic examples 

We digress to present a toy example of a discrete distribution of vevs in an iV-polygon 
enclosed by a ring of radius r in the 1-2 plane @ 

Xl ev = (r cos <f> p , r sin <f> p , 0,0, 0,0) , P = 2irp/N , p = 1, 2, . . . , N . (68) 

Computing the matrix elements using the definition (|67|) is straightforward. We find 
that the only non-zero components are An = A 2 2 = Ntq/2!^ We note that in this case 
we obtain the same result even if we approximate the discrete distibution by a continous 
uniform distribution of vev's on the circumference of the circle. 

We now turn to the specific examples of the distribution of vev's on a disc and on a 
three-sphere, which we considered in section 3 using AdS/CFT correspondence. In these 
cases a direct comparison with the free field calculation can be performed. In particular, 
in accordance with the convention in (|25|)-(|27|), the momentum space version of (^) can 
be expressed in terms of a function H(k) 

m ~ ' <™> 

For the distribution on a three sphere it is obvious that An = Nr^/A, for i = 1, 2, 3, 4 
and zero otherwise. These results are most easily derived in the continous approximation 
of the distributions. Hence, using (|70f) and the facts that g\ M = g s and i? 4 = 4iig s N, we 
obtain 

H(k) ~ -i * , (71) 

where the parameter A = 0, | and 1 corresponds to currents in the transverse direction 
(unbroken SO (2)), broken currents in the coset and directions along the distribution 
(unbroken £70(4)), respectively. This agrees nicely with the AdS/CFT result (|43|) . 

13 This is correct as long as we work with energies (distances in the gravity side) U not too close to 
the vev values. Typically the condition to be fulfilled for the continous approximation to be valied is 
U/X vcv - 1 > 0(l/N), where X vcv is a typical scalar vev value p|. 

1 We have used the fact that 

N N N 

cos 2 (2np/N) = Y sin 2 (2irp/N) = N/2 , ^ cos{2np/N) sin{2np/N) = . (69) 

p—1 p—1 p—1 
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For the distribution on a disc we have similarly that An = Ntq/4, for i — 1,2 and 
zero otherwise. Using ( fTO]) we compute 

r 2 \ — 1 

*(*) ~ ^ V , (72) 

where the parameter A = 0, \ and 1 corresponds to currents along the distribution (unbro- 
ken £0(2)), broken currents in the coset and directions orthogonal to the distribution 
(unbroken 5*0(4)), respectively. Again we find precise agreement with the AdS/CFT 
result (1581). 



Generalization to a class of models 

A natural question is whether the agreement between field theoretical results and those 
obtained from supergravity goes beyond the two specific examples we considered in detail. 
In fact, we may systematize our approach and show that such an agreement persists for all 
models with vev distributions corresponding to the five-dimensional supergravity solution 

(i)-(D- 

On the supergravity side the distribution of vev's is encoded in the harmonic function 
appearing in the ten-dimensional metric describing the gravitational field of D3-branes. 



In our cases the harmonic function is 11 



6 



where F is determined in terms of the six transverse coordinates yi as a solution of the 
algebraic equation 

Ef^ = 4 ' < 74 > 

i=i r Ul 

The harmonic function is in general 

N 



H = J2-^-, (75) 



where the vev values X p in (|65D became the centers of the harmonic function. In the 
continous approximation this takes the form 

H = Ang s [ (fx-^- , (76) 

J \y-x\ 

where the density p(x) is normalized as J d 6 xp(x) = N. We would like to compute A^ 
in (1671), which for a continuous distribution reads 



Aij = / d xp(x)xiXj . (77) 



22 



In general this can be found from the large r expansion 

» = fr - ^4 C*"- 5 ?^ (78) 

Returning to our cases where the harmonic function has the specific form (|73|) , we 
may cast its large r expansion into the above form with 

Aj = Nb Xj Sij , (79) 

where we define in general fey = bi — bj. We see that our general distributions allow 
a diagonal matrix Ay. Hence, the only non-zero independent components of the group 
theoretical factor H^* are H l ^ l K If all indices correspond to directions which are per- 
pendicular to the distribution then H^' = 0, whereas if all directions are along the 
distribution H^m = N(b\j + bu). If we are in the coset one index is along the distribu- 
tion (say i) and one is orthogonal to it (say j), then one of the above terms is missing 
and therefore H 1 ^^ = Nbu. This agrees perfectly with the two special cases of the disc 
and sphere distribution that we considered before. 



Correlators from supergravity 

Let us consider the equation (|2~3"D but in terms of the variable F 



d ( (F b)(F b) d ^ ^ (F-bi^F-bj) b% 
dF\ {F - bi){F - bj) dFj ~ k W 2 $ ~4(F-6,)(F-^ 



where F was defined in equation (g). Equation ( p0|) was solved exactly for the cases of 
the disc sphere and the sphere distribution. For the purposes of this addendum it suffices 
to concentrate on the limit k 2 — >• 0, where (]80f) can be solved exactly for any distribution. 
This will give the leading contribution to the two-point function of currents for large 
distances. At the AdS boundary F — » oo we impose the usual boundary condition 
$ — > 1 corresponding to a point-like source. Furthermore, we require $ to be smooth at 
the singularity F = b\ in the interior. In the following we use units where g = 2/R = 1. 

Currents transverse to the distribution : 

In this case the indices of the current i,j are such that 6j = bj = b\. Demanding 
regularity at the singularity F = b\ and imposing the normalization condition at the 
boundary gives 

$ = 1. (81) 



Therefore (24) gives 

H{k) = . (82) 
As expected this agrees with the field theoretical result. 
Currents longitudinal to the distribution: 



23 



In this case the indices of the current are such that 6j, bj =^b\. As before, regularity 
at the singularity at F — b\ and the normalization condition at the boundary give 

1 f 1 fF-bi\ 1/2 , fF-b^ 1 ' 2 



from which we compute using (pl[) that 



This is in perfect agreement with field theory expectations as spelled out after (|79|) . A 
particularly interesting case is when bi = bj ^ b\. Then the above expressions reduce to 

•=^£ < 85 > 

and 

H(k) = . (86) 

The case of the sphere and disc distributions correspond precisely to that with bu = Tq/4 
(b% can be put to zero by a shift of the coordinate F), for i = 1,2,3,4 and i = 1,2, 
respectively. 

Currents in the coset : 

In this case the currents indices are such that 6j = b% and bj ^ b\. Proceeding as 
before we find that 

and that 

m = -% ■ (88) 

Again, one easily sees that this agrees with field theoretical expectations. 



Comments on the masses of gauge bosons 

Finally, we mention some usefull facts about the masses of the W-bosons that arise on 
a generic point of the Coulomb branch of the J\f = 4 SYM theory. The general mass 



matrix is read off from eq. (63) of |46 



{M\ q =\X p -X q \\ p, q = 1, 2, . . . , N , (89) 

up to a numerical constant of order 1. Hence, the masses have the geometrical inter- 
pretation as the distances between the various vev positions distributed in the R 6 scalar 
space. Equivalently, they are given by the masses of the strings stretched between the 
D3-branes located at these points. Since we may shift uniformly all vectors X p 's without 
changing the Physics, the number of elements are in a generic case A^ 2 — 1 as it should 
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be. It is clear that depending on the specific vev distributions some of these masses 
might be degenerate. In particular, in the case of the discrete distribution of vev's in the 
iV-polygon we find, using fl68|), that (see eq. (66) of |46|| ) 



M n = r an{im/N) , n = 1, 2, . . . , JV , (90) 

which is an exact result for any N. The degeneracy for the zero mode is d,N = N — 1 and 
for the rest d n = 2{N — n). It is easily seen that Y2n=i — N 2 — 1. Hence, for large N 
there are W bosons with masses of order r and light masses of order r /N. In the case 
of vev's distributed on a disc a similar result can also be derived starting from a discrete 
distribution R [TT1 whose limit is the continuous one we have been using. 
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